Buckling instability in type-II superconductors with strong pinning 
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We predict a novel buckling instability in the critical state of thin type-II superconductors with 
strong pinning. This elastic instability appears in high perpendicular magnetic fields and may cause 
an almost periodic series of flux jumps visible in the magnetization curve. As an illustration we 
apply the obtained criteria to a long rectangular strip. 
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In high magnetic fields a noticeable deformation of su- 
perconductors occurs in the critical state because of the 
magnetic force density f = j x B, where j is the cur- 
rent density and B is the magnetic field. This results in 
an anomalous irreversible magnetostriction ( "suprastric- 
tion" 10]) and shape distortion j^^ of type-II supercon- 
ductors with strong pinning. Similar as in magnetic fluid 
dynamics [H the stress tensor of a superconductor in a 
magnetic field includes an additional term, the Maxwell 
stress tensor of the magnetic field with components of 
order B'^/^q. Since this is quadratic in B, the Maxwell 
stress tensor in the critical state can be important for 
elasticity in strong magnetic field y,|3|. However, even in 
a field of 10 T the value of B^ /^q is small compared to the 
Young modulus E of the material. We estimate the ratio 
B^/HoE « 10-3 for B « 10 T and £: « 100 GPa which is 
a typical value for YBa2Cu307_a; high-temperature su- 
perconductors g]. 

The effect of the magnetic field on the elastic behav- 
ior may be much higher if one considers bending of thin 
samples since the effective elastic modulus for bending E 
is much less than the Young modulus. In particular, for a 
long rectangular strip of extension lxwxd(l'^w^d) 
one has E « E{d/l)^ < E. If for instance d/l w 10"^ 
and E « 100 GPa, then B^/fio is of the order of the 
effective bending modulus E sA B k, 3.5 T. 

An important consequence of a small value of the effec- 
tive elastic modulus for bending E is the classical Euler 
buckhng instability |g,0. This elastic instability occurs 
for rods and thin strips when the longitudinal compres- 
sion force F at the edges of the sample exceeds a critical 
value Ff, ex E. In particular, one has F}, = ir'^Ewd'^/iSP 
for a long rectangular strip with one edge clamped and 
the other edge free as shown in Fig. |l| |^,Q. The buck- 
ling instability manifests itself at F > F;, by a sudden 
bending with amplitude s ex -^/F — Fj,. 

The magnetization of type-II superconductors with 
strong pinning and the associated magnetic forces are 
successfully described by the Bean critical state model 0] 



using a critical current density jc which decreases with 
increasing temperature and magnetic field. In the trans- 
verse geometry of a thin strip in a perpendicular field one 
has j = jc in the region where the magnetic flux has pen- 
etrated and screening sheet currents J with < J < jcd 
in the fiux-free region P-0| ■ This nonuniform flux distri- 
bution is not in equilibrium and under certain conditions 
a thcrmomagnctic flux-jump instability may occur pro- 
ducing a sudden intensive heat release. This heat pulse 
decreases the critical current density and drives the sys- 
tem towards the equilibrium state with a uniform flux. A 
sudden buckling of a superconductor in the critical state 
may also lead to a heat pulse and thus to a sudden flux 
penetration into the sample, which shows as a flux jump 
instability in the magnetization curve. 



F<F. 




FIG. 1. Buckling of a thin superconductor strip with a 
clamped left edge in a transverse magnetic field Ba ■ 

In this letter we predict a novel Euler buckling insta- 
bility caused by the longitudinal magnetic compression 
force acting in the critical state of a thin superconduct- 
ing strip in a strong transverse magnetic field. We discuss 
several scenarios how the buckling instability develops, 
including the cases when a sudden buckling shows as a 



flux jump instability in the magnetization curve. A se- 
ries of buckling induced flux jumps almost periodic in an 
increasing applied magnetic fleld is predicted. 




FIG. 2. The current streamlines in the critical state of 
a type-II superconductor thin strip in a transverse mag- 
netic field computed by the method [9]. The arrows indi- 
cate the magnetic forces acting on the strip. Top: Meiss- 
ner state, Ba <^ Be- Middle: Fully penetrated critical state, 
Ba^ Be. Bottom: Applied field decreasing from Bo ^ Be to 
-Bo — 2ABc, which yields penetrating fronts with inverse flux 
at |a;| — a/ cosh 2.4 ~ 0.56a and a negative force F, Eq. (0). 

We consider first the elastic stability of a long rect- 
angular strip / X w X d (^ 3> w S> d) in an increasing 
transverse magnetic field Bq Ijz, assuming that the strip is 
glued to the substrate at the left edge {y = 0) as shown in 
Fig. |l|. A longitudinal compression force F acts near the 
right edge of the strip {y — I) in the area where the elec- 
tromagnetic force density f = j x B has a y-component 
due to the U-turning current, thus 



F = -Fy = Bad / / jx dxdy, 



(1) 



where the integral is over the U-turn area. As shown in 
Fig. 0, in the fully penetrated critical state this area is 
a triangle where jy = jc, but in general the integral is 
over the right half of the strip {y > 1/2). If w ^ I the 
deformation of the strip can be obtained assuming that 
F is applied to the very end of the strip a,t y — I. For 
such narrow strips one can show that exactly F — BaM, 
where M is the total magnetic moment of the strip di- 
vided by its length I. 

Depending on the magnetic prehistory of the sample 



the dependence of the longitudinal compression force F 
on Ba is described by the following three formulas ||ll|,|l2[ . 
(a) For a zero-field cooled straight strip (Fig. la) with 
Ba increasing from zero one has 



F = jcBada^ tanh -^ , 

Or 



(2) 



where we introduce a = w/2 and Be = fJ-ojcd/T:. The 
longitudinal force F{Ba), Eq. (g), has the limits F « 
■KBaa'^/fJ.o {Ba ^ Be, Meissner state) and F w jcBada? 
{Ba ^ Be, fully penetrated critical state). 

(b) For Ba increasing from a field-cooled value Bq, one 
has the force 



F = jcBada^ tanh ^" ^° 

Br 



(3) 



(c) For Ba decreasing from a fully penetrated critical 
state with Ba — Bq, the force F = BaM decreases as 



F = jcBada^ 



1-2 tanh 



Bq ~ Ba 

2Br 



(4) 



going through F — at Ba f^ Bi^ — l.lBc. For a narrow 
strip the field of full penetration is fl^ 



Bp — Be 



1 w 

1 + ln- 

d 



(5) 



In the case of a curved strip (Fig. lb) in the formulae 
(2-4) for F the factor Ba = F/M means the z component 
of Ba, while in the argument of tanh(. . .) the Ba should 
be replaced by the component B±^ of Ba perpendicular to 
the strip near its right end (where the U-turning currents 
flow). In general, the magnetic moment M and the force 
F = BaM depend on the prehistory of B±{t) and may 
relax with time t. 

If the buckling instability for a zero-fleld cooled strip 
occurs at a field B^ > Bp, the force is 



F = jcda'Bb. 



(6) 



The critical force Fi, for the buckling instability of a_s_trip 
with one edge clamped and the other edge free is 



Fh = 



24/2 



(7) 



Equating the forces F and Ff, we find that the magnetic 
field Bb at which the first buckling instability occurs is 



Bh^ 



Fb 
jcda^ 



TT^ E /dy 
24 iTa Uy 



(8) 



We estimate the fields B^ w 0.04 T, Bp « 0.15 T, and 
Sft w 4T using the data for YBa2Cu307_a; superconduc- 
tors E w 10^ GPa |||, and assuming that jc ~ 10^ A/m^, 
w « 10"^ m, d « 10"" m, and d/l « 10"^ This estimate 
verifies our initial suggestion that Bp <^ Bb. 



The height s of the right end of the buckled strip (see 
Fig. ^3) can be found analyticaUy if the maximum angle 
9m between the tangent to the buckled strip and the sub- 
strate is small ||,0|. Assuming that the force F slightly 
exceeds the critical value Fb we obtain a sinusoidal bend- 
ing with the amplitude 



4V2 



F 



'■^4 



and 



I- I F TT S 



(9) 



(10) 



Now assume that the external magnetic field is in- 
creased with constant ramp rate Ba and the threshold 
of the buckling instability is reached when F — Fb- One 
can consider several scenarios how the buckling evolves, 
depending on the value of Ba and on the ratio of the time 
constants for bending of the strip, Tb, for magnetic flux 
diffusion, r,„, and for heat diffusion, t/j, see Ref. |13| for 
details. 

The first scenario applies to a very low ramp rate 
Ba <C BaJTm, whcrc the current and magnetic field dis- 
tributions inside the strip, and thus the magnetic forces, 
follow the increasing field Ba without delay. In this 
case the strip starts to bend as soon as the magnetic 
compression force F reaches the critical value Fb- The 
force F = BaM via the magnetic moment M depends 
on the perpendicular field component near the tilted tip 
of the long strip, B±^ = Ba cos 9„i- This means that in 
9m{F), Eq. (p^), F depends on 6',„ and one has to find 
the value of 6jn self-consistently. To do this we need 
the appropriate dependence M{B±^). We shall see that 
the resulting B± decreases with increasing Ba (or time); 
thus we have to use Eq. (Q) with Bq = Bb (the field 
where buckling starts) and with Ba replaced by B± in 
tanh(...). Expanding the hyperbolic tangent we thus 
find for Bb-Ba<.Bc, 



F^Fb 



Bb 



1 



Bb — B± 

b7~ 



(11) 



Inserting this force into Eq. (p^), 0^„ = 8{F/Fb — 1), and 



solving for 
we obtain 



using B± « Ba{l - 61^/2) and Ba > Be 



V2 



B^ 
Bb 



1. 



(12) 



This self-consistent tilt angle Om is two times less than 
the tilt angle Eq. (HO) for constant compression force F. 
The physical origin of this negative feedback is the reduc- 
tion of the total U-turning current and thus of the force 
F, caused by the decrease of B±^ when the end of the 
strip tilts, compare the current distributions in Fig. El 



A different scenario appears when the buckling occurs 
with a delay at a force Fd slightly above Fb ("overheat- 
ing"). Several reasons for such a delay are conceivable, 
e.g., sticking of the strip to the substrate by adhesion, 
or a misalignment of the perpendicular applied magnetic 
field Ba such that the force F in Fig. ^ points slightly 
downward to the substrate. A small misalignment is 
probably inevitable for a typical experiment. 

When after zero-field cooling F = Fd = jcBdda^ is 



reached at Ba 

most instantly to a finite value s 



Bd, the buckling amplitude jumps al- 



1. To 

obtain this amplitude self-consistently one may combine 
Eq. (0) for 6'™(F) with Eq. (|) for F{9m), like in the first 
scenario, noting that M and thus the force F = BaM 
depend on B± = Ba cos 9^- The sudden jump of 9^ 
at Ba = Bd means that B± is reduced from Bd to 
Bd{l - 91^/2) (if 91^ < 1) and thus Eq. (|) is required 
yielding 



F^Fd 



1-2 tanh ■ 



ABr 



(13) 



and 



1 



(14) 



with Fd = jcBdda?. Inserting this into Eq. (|1C 
solving for 6'^„ <C ABc/Bd one obtains for Bd ^ Be 

2 ^ 2_Be / Fd 
™ Bd VFfc 

Equation (O) differs from Eq. (12) because of the dif- 
ferent history of the perpendicular field B± (t) and thus 
of the magnetic moment: In the first scenario B± started 
to decrease from the lower threshold Bb and the decrease 
occurs since the rising Ba is overcompensated by the 
growing 9m- In the present scenario, B± has reached 
the higher threshold Bd > Bb before it drops down, and 
this drop is solely due to the growing tilt angle 6',„ while 
Ba ~ Bd is constant in this approximation. As a conse- 
quence, the self-consistent tilt angle 9^ is reduced much 
more in this case, by a factor ^jBc/4,Bd ^ 1. 

This strong feedback mechanism requires that the 
change of the current density occurs instantaneously^ 
much faster than the mechanical buckling, t,„ ^ r^. 
In reality the redistribution of the currents will lag be- 
hind the buckling. In the extreme limit t,„ ^ r^, the 
tilt angle would first jump to its original large value 
91 = 8{Fd/Fb - 1), Eq. ^, and then relax to the 
small value of Eq. ([l4|), or to zero, or to some other 
value. The theoretical problem is intricate since a quan- 
titative treatment requires the self-consistent time de- 
pendent solution of the equations for Bj_{t) with a re- 
laxing, history dependent magnetic moment M{B±{t)}, 
using Bx = Bait){l - 91/2) and 9^ = 8iF/Fb - 1) with 
F — BaM. This yields the implicit equation for B±{t), 

B^{t) = Ba{t)[5 - 4Ba{t)M{B^{t)}/Fb] , (15) 

from which the tilt angle 9l^{t) ^ 2{B^/Ba - 1) is ob- 
tained. To solve this one requires a realistic model for 
the relaxing history dependent magnetization. 



From our numerical work we expect the magnetic 
relaxation to be very fast and non-exponential when 
dB±/dt changes sign Q, as it is the case during buck- 
ling. During very fast switching of B±{t), the electric 
field is so large that, irrespective of pinning, the vortices 
exhibit usual flux-flow behavior, with flux-flow resistiv- 
ity pf ~ {B/Bc2)pn, where Bc2 is the upper critical field 
and pn is the resistivity in the normal state. In this 
case the magnetic relaxation time of an Ohmic strip ap- 
plies, Tm ~ To = 0.249adpo/pf jl^]. This time has to 
be compared with the buckling time rb, which we esti- 
mate from the lowest resonance frequency ui of the strip 
(a cantilevered reed H), n ~ uji\ ujf « 1.03 Ed^/{pl^) 
where p is the specific weight. Inserting here numbers for 
YBa2Cu307_x a.t Ba — 4: T, we estimate r™ ^ Tb, i.e., 
the magnetic relaxation initially is instantaneous. With 
proceeding relaxation, the electric field and the effec- 
tive resistivity decrease, and thus the magnetic relaxation 
time increases. We thus expect that the real behavior of 
the strip is somewhere between the two considered limits 
r„i < n and Tm > n. 

Therefore, if buckling starts delayed at a force Fd and 
disappears at a smaller force Fb < Fd, the tilt angle at the 
tip of the strip may oscillate between a maximum value 
^max < dd and zero. Such oscillation may occur since at 
dd the reduction of B± is so large that the currents tend 
to change sign and thus the force F rapidly decreases. 
The tilt angle then may drop to zero, undershooting the 
small equilibrium value, Eq. (p^). With continuously in- 
creasing applied field Ba{t), the tilt angle Om thus makes 
a sudden jump from zero to ^max, then drops rapidly 
back to zero, where it remains until the next excursion 
occurs when F again reaches Fd. These buckling instabil- 
ities should occur at nearly equidistant field values with 
period of the order of the penetration field Bp, Eq. (||), 
and they will show up in the magnetization curve as a 
periodic set of fiux jumps. 

So far we assumed that the temperature T of the strip 
stays constant, T ^ Tq. However, buckling of a strip in 
the critical state causes some heat release which increases 
the temperature and decreases the force F{T) ex jc{T). 
A complete solution of the buckling instability in type-II 
superconductors with high critical current density should 
therefore include a self-consistent treatment of the mag- 
netic field and temperature variations. 

For a rough estimate of the decrease of the force F{T) 
we assume here that jc{T) ex (Tc — T), the critical tem- 
perature is Tc^ To, and the heating of the strip is adia- 
batic. In this case wc find that a sudden tilt to an angle 
6m leads to 



F{T) ~ F{To) 



jcwBd 



C{To)Tc 2 



(16) 



C{Tq)Tc <S^ jcwBc, which results in Tq ^ 3K for a heat 
capacity C{T) « 7 x T^J/Km^ ||]. 

The temperature dependence of F{T) may cause oscil- 
lations of the strip. Indeed, a sudden buckling leads to a 
heat pulse increasing the temperature T and decreasing 
the force F{T). If because of the temperature increase 
the force F{T) falls below the buckling threshold Ff, then 
the strip straightens and the next instability occurs after 
the strip has cooled down. 

In summary, we have shown that a strong magnetic 
field applied perpendicular to a cantilevered supercon- 
ductor strip will lead to Euler buckhng of this strip. We 
give the threshold field at which this clastic instability 
occurs. During buckling, the effective applied field at 
the tip of the strip decreases due to tilting. As a con- 
sequence, the buckling force is reduced. This feedback 
mechanism may lead to mechanical oscillations of the 
strip and its magnetization, which depend on the mag- 
netic and thermal relaxation times of the specific experi- 
ment. At sufficiently low temperatures this sudden buck- 
ling may trigger a periodic series of flux-jump instabilities 
which should show in the magnetization curve. 
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Combining Eqs. ®, (|l3|), and (|l|) we find that 
self-heating affects the buckling instability threshold if 
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